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Abstract. We investigate the properties of two- and three-dimensional non- 
commutative fermion gases with fixed total z-component of angular momentum, Jz, 
and at high density for the simplest form of non-commutativity involving constant 
spatial commutators. Analytic expressions for the entropy and pressure are found. 
The entropy exhibits non-extensive behaviour while the pressure reveals the presence of 
incompressibility in two, but not in three dimensions. Remarkably, for two-dimensional 
systems close to the incompressible density, the entropy is proportional to the square 
root of the system size, i.e., for such systems the number of microscopic degrees of 
freedom is determined by the circumference, rather than the area (size) of the system. 
The absence of incompressibility in three dimensions, and subsequently also the absence 
of a scaling law for the entropy analogous to the one found in two dimensions, 
is attributed to the form of the non-commutativity used here, the breaking of the 
rotational symmetry it implies and the subsequent constraint on J^, rather than the 
angular momentum J. Restoring the rotational symmetry while constraining the total 
angular momentum J seems to be crucial for incompressibility in three dimensions. 
We briefly discuss ways in which this may be done and point out possible obstacles. 



PACS numbers: ll.lO.Nx, 05.30.Fk, 02.40. Gh, 05.70.Ce 



1. Introduction 



The realization that gravitational instabilities make the absolute localization of events 
impossible, which in turn may be captured through a non- commutative space-time 
structure p], has led to considerable interest in the formulation of non-commutative 
quantum field theories |2] , quantum mechanics [3] and the possible physical consequences 
of non-commutativity |2|5l|6l[7l[8l|9l[l0l[T3l[l2l[Tl]. 

One aspect of non-commutative systems that has received very little attention 
in the literature is the thermodynamics. The reason for this is probably three-fold. 
Firstly, to discuss the thermodynamics one has to consider systems confined to a 
finite volume, e.g. a box, to investigate the dependence of the thermodynamics on 
the system size. The introduction of boundaries for non-commutative systems is 
notoriously difficult and has only been done successfully for a two-dimensional well 
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[IS] . Secondly, apart from the example just mentioned, there are no other non-trivial 
(apart from the harmonic oscillator) systems for which the energy spectrum is known, 
even numerically, which makes the computation of the thermodynamics intractable. 
Here it should be mentioned that the spectra of a number of non- commutative systems 
have been computed perturbatively to leading order in the non-commutative parameter 
UlElEllTllHlElITOllISlllSlIll]. However, as will become clear later, this is not sufficient 
if one wants to understand the physically interesting high temperature or high density 
limits, which show strong deviations from commutative behaviour. Thirdly, the issue of 
the breaking of rotational symmetry in three dimensions creates a conceptual barrier. 
One either has to accept it, seek modified commutation relations that respect the 
symmetry or one can restore the symmetry through twisting [131 [IHl EZl llHl ttS]. If 
the latter route is chosen, it also becomes necessary to twist the statistics, but there 
seems to be disagreement on whether this should be done in the sense of [HI [19] , which 
results in modified statistical correlations [19j , or in the braided sense of [201 121] • 

The first attempt to understand the thermodynamics of such systems was made in 
|22j . In this case the difficulties above were avoided by considering particles confined 
in a two-dimensional infinite well for which the spectrum is known [15]. The issues of 
rotational invariance and the twisting of statistics are also not present. However, the 
computation could still only be carried out numerically and thus the thermodynamic 
limit could not be probed reliably. The results of these calculations were quite 
striking and showed strong deviations from commutative behaviour at high densities 
and temperatures. The most important feature is the appearance of incompressibility 
at high densities for fixed angular momentum. This refiects the effective excluded 
volume implied by the non-commutative length scale. Under these conditions several 
other thermodynamic features also deviate strongly from the commutative case and, in 
particular, the entropy is no longer extensive. 

The surprising nature of these results provide a strong motivation for a better 
understanding, both in terms of a more reliable analytic treatment that can probe the 
thermodynamic and high density limits more efficiently and the generalization to three 
dimensions. With regard to the latter, the issue of symmetry breaking still prevails. Here 
we simply accept the breaking of rotational symmetry and investigate the consequences 
without any attempt to restore the symmetry through modified commutation relations 
or twisting. The naive expectation is that the breaking of rotational symmetry should 
be irrelevant in the thermodynamic limit. The essential reasoning behind this is that 
the shape of the container we use to study the system, be it a sphere or a cylinder, 
is immaterial in the thermodynamic limit and only shows up in finite size corrections. 
However, it turns out that this expectation is only realized in the low density limit, while 
it fails in the high density limit. The rest of the high density thermodynamics also do 
not reflect the desired features and in particular incompressibility. The expectation, 
therefore, is that symmetry breaking matters at high density and that the restoration 
thereof will be a key element for the understanding of the high density thermodynamics. 
The possible scenarios for restoration of the rotational symmetry and the possible 
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problems associated with them are discussed later. 

This sets the scene for the present paper, which aims at understanding the 
thermodynamics of non-commutative fermion gases in two and three dimensions 
analytically, particularly at high densities and low temperature. The paper is organized 
as follows. Section [2] reviews and generalizes the results of |22] for the two-dimensional 
gas. The focus here is on the analytic treatment of the high density limit. Section 
12.51 focuses on the two-dimensional thermodynamics close to the incompressible point 
and it is shown that the entropy exhibits a square root dependence on the system size. 
Section [3] generalizes these results to three dimensions and discusses the issues around 
the restoration of the rotational symmetry. The paper is concluded in section H] with a 
summary and outlook. 



2. The two-dimensional non-commutative fermion gas 

2.1. Definitions 

We begin by summarizing the key results of [IHl [22] . Consider a system of spinless 
particles of mass tuq which are confined to a disc with area A = ttB?. The particle coor- 
dinates satisfy [x, y] = iO {0 > 0) and we set = 6{2M + 1) with M a positive integer. 
The latter serves as a dimensionless measure of the system size. The single particle 
energies are then Er^m = f^Xr^m/ {d^o) with Xr,m the zeros of the Laguerre polynomial 
-^M+i(^) ("^ — 0) -^M+m+i(^) < m < 0). Here m is an angular momentum 

label while r labels the M + 1 (m>0)orM + m + l (— M < m < 0) states within a 
specific angular momentum sector. Note that m is bounded from below by —M and so 
there are only a finite number of negative angular momentum states. 



The g-potential of the grand canonical ensemble is 

oo 

q{A, P,fi,u)= J2 Yl + e-^^'^^'-^-f'-^"'^] . (1) 

m=—M r 

We set Eq = h^/{9mQ) and introduce the dimensionless parameters (3 = Eof3, jl = h/Eq 
and uj = fvjj / Eq in terms of which the g-potential becomes 

oo 

g(M,/3,/i,(:;)= log[l + e-'^^^-"^"^'")]. (2) 

m=—M r 

The central thermodynamic quantities can be calculated from q using 

.r L Idq S ~dq 

h pdoo k dp 

and we also define dimensionless measures of the density and pressure by p = N/ (2M+1) 
and P = g/(/3(2M + 1)). 



The entropy of dense nan- commutative fermion gases 







\f(Q')=c+wa 






6 










2 








-1 


a. 1 


1 1 1 1 1 TV 1 1 1 

2 "+ \ 3 





/ 2.5 






le fry) 9 






/ 1.5 




S 


j 1.0 


✓ 

✓ 


fi£)__JN 


Y 0.5 




-2 


-1 


\ «+l 2 



(a) 



(b) 



Figure 1. The curves e±(a) define tlie range of tire rescaled single particle energy e for 
each value of a. These two curves meet at the point (—1, 1) which is indicated by a dot. 
(a) The line /(a), shown here for c = 7.4 and w = —3, defines an effective chemical 
potential within each angular momentum sector. The points a_ and a+ denote the 
intersections of /(a) with the e±(a) curves. At low temperatures the integrand in (|6]) 
is appreciable only within the shaded region, (b) An illustration of the low density 
case with small values of c and w. 



In the thermodynamic hmit ([2]) can be replaced by 

/oo px^lm) 
dm / da; D{x, m) log[l + e" 
-M J x^{m) 



(4) 



where D{x^ m) is the asymptotic density of the zeros of the Laguerre polynomial L^^^(a;) 
(m > 0) or l}^'^j^^_^i{x) (— M < m < 0) in the limit where M and m tend to infinity in 
a fixed ratio. It is known that D{x,m) is given by [21] 



D(x, m) 



iy/4Mx — (m — xf 
2ttx 



(5) 



for X in the interval on which the argument of the radical is non-negative. 
2.2. The low-density limit 

To analyse the integral in (j4j) we first transform to the rescaled parameters e 
a = m/M, b = M/3, c = jl/M and w = u. The g-potential then reads 

r.e+(a) 



da / de rf(e, a) log[l + e" 

1 Je-{a) 

where /(a) = c + wa and 

^47 



-6(e-/(Q))i 



d(e,a) = D{Me,Ma) 



[a - e 



27re 

The latter is defined for e G [e_(a;), e+(a)] with 



e^(a) = a + 2±2V« + l. 



(6) 

(7) 
(8) 
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The form of the integrand in suggests that /(a) acts as an effective chemical 
potential within the m = aM angular momentum sector. At low temperatures al- 
most the entire contribution to the integral therefore comes from the region defined by 
< e < min[/(a), e+(a)]. This is illustrated in figure [1] for particular choices of 
parameters. 

Next we illustrate how at low densities non-commutative effects become 
unimportant and the physics reduces to that of a standard commutative fermion gas. 

Consider a process in which the system size M is increased at a fixed particle density 
and angular momentum L = 0. Since the chemical potential jl is intensive we expect 
c = fi/M to approach zero as M tends towards infinity. From figured] (b) it is clear that 
to produce a zero total angular momentum w must itself be approximately zero with the 
intersections a_ and a+ located roughly symmetrically around a = at a± ~ ±2^^. 
At low temperatures almost the entire contribution to ([6]) will therefore come from a 
region where |a| ~ y/c and e ~ c. This implies that in (E]) we may approximate d{e, a) 
by 



a 



2 



d{e,a)^ (9) 

where only terms up to order ~ has been kept in the square root. This is 

equivalent to approximating the density of zeros in ([5]) by D{x, m) ^ \/4Mx — w? / {2'kx) 
or the density of states in the m angular momentum sector by 



D(E) « V^^'"°^;ff - (10) 

It remains to show that this is precisely the density of states of a commutative 
fermion gas confined to a disk with radius R. For such a gas the single particle energies 
are given by i?j.,m = ^'^jrm/i'^^oR'^) where jr,m are the zeros of the Bessel function 
Jm{jr,m) = with 771 E Z and r G N. The zeros {jr,m} are distributed with density [23] 

D,ov-«) = ^iH^ (11) 

and it is straightforward to show that this indeed produces the density of states in ( ITOll . 



2.3. The high-density limit 

The results of the previous section suggest that to observe non-commutative effects in 
the thermodynamic limit we must consider densities which scale with the system size. 
This will render the chemical potential extensive and ensure that c = fl/M remains 
finite as M tends to infinity. A more appropriate measure of the density is therefore 
u = p/M ^ N/{2M'^), which will be kept fixed when taking M — )■ oo. We also define 
the scaled angular momentum by £ = L/{M^h). 
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Figure 2. Shaded regions indicate the physically allowed values of the scaled density 
v — p/M and angular momentum i = L/(M^h). Subfigure (b) is an enlargement of 
(a) around the origin. 



The integrals defining the g-potential can be approximated at low temperatures 
using the Sommerfeld expansion |25]. Appendix A contains a summary of this 

: order in temperature, reads 



procedure. The final result, up to linear order in temperature, reads 



M2 

6b 



[w 



q{M,b,c,w) / 

where 0(-) is the step function. This result is valid for w < 1 and c > w"^/ {w — 1). The 
first condition is necessary to ensure that the integrals in ([6]) converge. If the second 
condition is violated then the /(a) line does not intersect the e±(a) curves (see figure [I]) 
and so g(M, b, c,w) = within this approximation. 

Transforming from {b,c,w) back to {P,fl,u}) yields the desired expression for 
q{M, 13, fi,uj) in ([2]) in the high density, low temperature limit. 



2.4- Solving for ji and u 

Our next task is to combine ([3]) and ( |T2i) and solve for fl and cD as functions of the 
angular momentum and number of particles. It is convenient to consider the three cases 
(I) A > 0, (II) ft < 0, oj > and (III) fl < 0, u < separately. Each of these corresponds 
to a physically allowed region of the {u, £)-plane appearing in figure [2l The three regions 
are bounded by the curves 

Ci(zy) = 8z/^/V3 - 2z/ C2{iy) = v{2v - 1) - 1/24 (13a) 

C;{u) = 8z/2 - I6Z/3/73 C^{u) = 16z/=^/V3 + 8z/2 (136) 
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with Ci{u) and Cs{u) defined for < z/ < 1/4, Ciiu) for z/ > 1/4 and C4(z/) for u > 0. 
We see that the angular momentum is bounded from below by —M^/6; a reflection 
of the finite number of negative angular momentum states. A more striking feature 
apparent in figure [2] is that for each i > —1/6 there exists a maximum density pc = MUf. 
with 

z/c = (3 + 2Vl8£ + 3)/12. (14) 

No such restriction is present in the commutative case and this is a clear indication that 
non-commutativity results in an excluded area effect at high densities. We investigate 
the thermodynamic properties close to this critical density in the next section. 



In regions (II) and (III) exact solutions for p and u can be found in the 
thermodynamic limit. The results are 



M 



V 



3£ + 4 (3 ± 507) V 



6u 



nl/2 



Se ± 8Z/3/2 ^ 6 



(15) 



[6z/(3£±8z/3/2 + 6z/)]'/' 

where the top and bottom signs correspond to regions (II) and (III) respectively. In 
region (I) the exact expressions for p and u are quite cumbersome and so we will focus 
only on the physically interesting region close to the critical density. Here we find 



A = l 
M 



7 + 24£-24z/(6z/- 1) 
24 [(8z/ + 2)(z/-z/o)(z/e-z/)]^/=^ 



1 

1-4 



(z/ - z/o)(z/c - z/) 



1/2 



(16) 



where z/q = (3 - 2^/l8iT3) /12. Note that as z/ approaches Uc these quantities diverge: 
/i — )■ oo and u — )■ — oo. In figured] (a) this implies that the "Fermi-line" /(a) = c + wa 
becomes vertical at the critical density and that all states in angular momentum sectors 
to the left of it are occupied. 



2.5. Thermodynamics close to Vc 



Substituting the expressions for p and u in ( IT6l) into the relations ([3]) yields expressions 
for the entropy and pressure of the non-commutative gas close to the critical density: 



S 



l^c - 1^0 

l + 4z/c 



1/2 23/2 



TT 



3/3 



and 



P 



l + 4z/, 

- 



.1/2 



1 + 48z/2 



96^2 V'^c - V 



,(17) 



In terms of the dimensionless system size M and particle density p ~ N/ (2M) this 
implies, for example, that at £ = 



s ~ Vm^/m^c 



p and P~ M^'^ /^Mvc- P- li»J 

We see that as z/ approaches Uc the entropy vanishes as the square root of Uc — z/, 
indicating that the "tight packing" limit has been reached and only a single microstate 
is accessible to the system. This coincides with a divergence in the pressure which 
signals that the system is incompressible at the critical density. Due to their entropic 
nature these phenomenon are expected to persist at higher temperatures. In particular, 
the value of Uc and the i/z/c — z/ dependence of the entropy and pressure should remain 
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(c) (d) 

Figure 3. Solid and dashed lines correspond to non-commutative and commutative 
results respectively. For the entropy plots f3 is arbitrary but should be sufficiently 
large {(3 > 100). (a) The entropy versus at ^ = 0. In the commutative case 
S/k = M7r^/(3/3). (b) The entropy versus M for TV = 10*^ particles and ^ = 0. Here 
Mc is such that i^c — N/{2M^). (c) The entropy versus at ^ = 1.2. (d) The pressure 
versus v a,t £ — for f3 — 100. The temperature dependence of P is however very weak. 

unaffected provided that the system is sufficiently close to its maximum density. Figure 
|3] (a) and (c) show the entropy as a function of u ai i = and i = 1.2. The cusp in the 
entropy in the second case occurs at the crossing of the C4 curve from region II to I. 
(See figure O) In figure |3](b) the number of particles is kept fixed and the system size 
increased from its minimum value of M^. Here we again see a clear deviation from the 
purely extensive behaviour of the commutative case. Figure [3] (d) shows the pressure as 
a function of u with the square-root singularity appearing at the critical density. 

The next question we address is the following: How will the entropy of systems 
close to the incompressible density Uc scale with system size? The physical motivation 
for this question is that we would like to know how the entropy of a very dense gas, close 
to its maximum density, behaves as a function of the size or number of particles when 
more particles are added and the size of the object increased, but in such a way that 
the density is always kept very close to the maximal density. For this purpose consider 
a system containing N particles with a minimum size given by = N/{2i>c). The 
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density can then be expressed as u = N/{2M'^) = UcM^/M^. The hmit that interest 
us is when the system size exceeds the minimum size by a small fixed area 27c6AM. 
Then M = Mc + AM with AM a non-negative integer such that AM << M^ and 
p ~ z/c(l — 2AM/Mc), i.e., the density is close to the incompressible density. Crucial 
here is that the quantization of the system size implies that AM ~ 1. With AM held 
fixed we now consider the dependence of the entropy on the system size M M^.. 
Substituting in (fT7|) the entropy close to the incompressible point reads 



As B? ^ 2M9 it follows that the entropy is proportional to the circumference, rather 
than the area, of the system. When M » Mc, i.e. far from the incompressible point, 
normal extensive behaviour is recovered. 

Interestingly, the proportionality constant takes discrete values (AM = 0, 1, 2 . . .). 
The way to interpret this is probably to think in terms of families of systems 
parameterized by the number of particles or the minimal system size. Each family 
is characterized by a specific value of AM, i.e., for minimal size AM = the entropy 
vanishes, for AM = 1 the entropy behaves as S{Mc) ~ a/2z^a/Mc etc. Within the 
current context it is, of course, not possible to determine the value of this proportionality 
constant. 

This is a rather remarkable result, especially since it follows from a first principle 
calculation, as it implies that the number of microscopic degrees of freedom for such 
dense systems is actually determined by the circumference rather than the area (size) 
of the system. It is also worthwhile noting that the argument that leads to this result 
has only a number of generic features and essentially relies only on the existence of 
an incompressible point and the quantization of system size. This robustness of this 
result can also be understood from an intuitive point of view. In essence the minimum 
size corresponds to "tight packing" for which only one microstate is possible and thus 
the entropy vanishes (see figure 1). For values of the system size slightly above than 
the minimal size, and thus for densities slightly lower than the incompressible density, 
only a fraction of the particles in a thin shell around the Fermi level determine the 
thermodynamic behaviour; the bulk of particles is still frozen in a "tightly packed core" . 
Furthermore, the particles in this thin shell occupy the highest single particle angular 
momentum states and are therefore expected to be close to the edge of the well. From 
this the intuitive picture is therefore that in real space a thin shell of particles at the edge 
of the well determines the thermodynamic properties of the system. Indeed, assuming a 
shell thickness AR and that each particle occupies an area 9, one easily concludes that 
the number of microstates behave as i7 ~ 2'^'^^^^/^ and thus that the entropy behaves 
as S ^ R with the proportionality constant set by the shell thickness as was also found 
above. From this perspective the behaviour found above is not surprising and expected 
to be rather generic. 




(19) 
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3. The three-dimensional non-commutative fermion gas 

We assume the following commutation relations in three dimensions 

[xi,Xj] = i9ij (20) 

with 6ij a real anti-symmetric matrix. Through an appropriate choice of coordinates it 
is always possible to restrict the non-commutativity to two of the spacial coordinates 
with the third coordinate being commutative. We consider a cylindrical geometry with 
radius R and length L with its central axis orientated along the z-direction. The x- and 
y-directions are taken to be non-commutative and at a fixed z the problem therefore 
reduces to the two-dimensional case considered earlier. In the z-direction we have the 
regular one- dimensional free particle problem with vanishing boundary conditions. In 
the notation of ([2]) the g-potential reads 

oo 

g3D(M, P, fi,u,^)=Y,Yl + e-'^(^'--+^"'-^-^'-)] (21) 

n=l m,r 
oo 

= J2<lMMj,t^-W,Cj) (22) 

n=l 

where 7 = 7/-E0 and 7 = h'^Ti'^ / {2moL'^) . In the thermodynamic limit R and L will 
tend to infinity in a fixed ratio. This can be made explicit by writing ^ = G/M where 
G = (7r/2)^(i?/L)^ is a dimensionless constant defining the shape of the cylinder. As 
is evident from (122|) the gs^j-potential is just a sum of two-dimensional potentials with 
shifted chemical potentials. The sum terminates automatically once jl — ^n^ drops below 
MCj'^ /{Cj — 1), as discussed after equation f[T2|) . Under the same assumptions that led 
to q2D in ( IT2I) we can now express q^n as 

q3D{M,P,fi,u,^) = dn q2D{M,(3,p, -^n'^,Cj) (23) 
Jo 

where n+ > is such that /i — 7^^ = MCd"^ /{Co — 1). This integral is straightforward to 
perform since g2D is a piecewise polynomial in its third argument. 

3.1. Scaling behaviour of the entropy 

In three dimensions it is convenient to use v = N/M^ and i = L/ (HM^) as measures of 
the density and angular momentum. Combining f l23|) and (|3]) leads to a set of equations 
for ft and u which can be solved numerically for given u and i. This produces the 
results of figure HI In (a) we show the entropy as a function of i' for i = 0. Unlike the 
two-dimensional case there is no maximum density at which the entropy vanishes and 
pressure diverges, i.e. no incompressible point. This refiects the fact that the excluded 
area effect is restricted to the x — y plane and does not place any restriction on how 
the commutative modes in the z-direction are occupied. However, the effects of non- 
commutativity are still present and responsible for modifying the scaling behaviour of 
the entropy as we progress from low to high densities. This is clearly reflected in flgure 
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m (b) which reveals that S grows hke u^^^ at low densities and decreases like 1/z/ at high 
densities. In these two limits we are able to solve for fi and u when £ = to find 

fi = M (9GZ/V4) u = - (18GZ/2) /5 (^^ « 1) (24a) 

fl = M {0.72703... )Gu^ u = -{0.823A2 . . .)Gu^ (u » 1) (246) 

The constants in fl24&p are derived from the numerical solution of a transcendental 
equation. The entropy is 



S 



TT 



6 ni/3 



1/3 ^ S (4.2397...) 1 
~ J' and TiTJ^ = ^- (25) 

for low and high densities respectively. These results provide the dashed lines in the 
log- log plot in figure |l](b). In terms of the dimensionless volume V = V/O^^"^ and density 
p = N/V these expressions become 

- = - l^pV3 and - = -!-^ . 26 

with 7] fa 6.4757 X 10"'^. The low density result agrees exactly with that of a commutative 
fermion gas [25j. It is independent of both the shape of the cylinder (i.e. G) as well 
as the non-commutative parameter 6. In contrast, the high density result depends on 
both these parameters and also exhibits an unusual non-extensive volume dependence. 
Furthermore, this entropy is inversely proportional to the particle density as is the 
case for a one-dimensional ideal gas. This observation that some aspects of the high 
density thermodynamics mimic that of a one- dimensional system is also observed in the 
relation between pressure and energy density. For a non-interacting gas in d dimensions 
it holds quite generally that P = {2/d){E/V). We indeed find that P = {2/3){E/V) 
and P = 2{E/V) in the low and high densities limits respectively. 

These results hint that this type of non-commutativity and constraint on lead 
to a high density limit that behaves as a commutative one-dimensional system. This 
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counters our expectations for the real world and we can therefore only assume that this 
form of non-commutativity, the symmetry breaking it implies and the constraint on 
are inappropriate. In particular the dependence on the shape of the container, i.e., the 
G dependence above is apparently due to the breaking of rotational symmetry. This 
is confirmed by comparing with the commutative case where such a dependence also 
occurs due to the constraint we imposed here, i.e. due to the breaking of rotational 
symmetry. This suggests that the restoration of rotational symmetry and the imposition 
of the more appropriate constraint on J are crucial for recovering the appropriate high 
density behaviour. In this regard it is important to remark that we do not expect 
that the restoration of the symmetry through twisting and the corresponding issue of 
twisted statistics alone can lead to the appropriate high density thermodynamics, but 
that the constraint on J is also a key ingredient. The reason for this is that the twist 
simply implies a unitary transformation implemented through the twist operator \18\ 
and thus should not affect the thermodynamics as the computation of the partition 
function involves a trace. However, correlation functions are affected by twisting as was 
shown in [T9 j. 

Attempts at restoring the rotational symmetry through twisting and constraining J 
may also be problematic as it was shown in [281129] that even though twisting restores the 
vectorial transformation properties of the non-commutative coordinates, the symmetry 
is still not manifest on the level of the Hamiltonian and Schroedinger equation, i.e., the 
Hamiltonian does not commute with J^. It is therefore not clear that this constraint 
can be imposed consistently in a statistical ensemble setting. 

A natural alternative is to change the non-commutative relations from the simple 
form (j20l) to those of a fuzzy sphere = iOeijkXk- In this way can be realized 

as an "onion structure" with quantized radius and the two-dimensional program of 
[22] can in principle be executed to compute the spectrum of a three-dimensional well. 
Unfortunately this program is technically extremely demanding and has not yet been 
performed successfully. Some progress in this direction was made in [26] and seems to 
confirm the expectations described here. 

4. Conclusions 

The thermodynamic behaviour of two- and three-dimensional non-commutative gases 
was investigated analytically in the high density region. Strong deviations from 
conventional commutative behaviour was found. In particular an incompressible point 
was found in two dimensions and the entropy exhibits in both cases non-extensive 
behaviour. More profoundly, for two-dimensional systems with density close to the 
incompressible density the entropy exhibits a very particular dependence on the system 
size in that it is proportional to the square root of the system size, i.e. the circumference 
rather than the area. Incompressibility is absent in three dimensions and the scaling of 
entropy at high density is unconventional, showing rather a crossover to one- dimensional 
commutative thermodynamics. This is probably due to the particular form of non- 
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commutativity used here in three dimensions and the symmetry breaking it imphes. In 
particular the constraint on J^, rather than J, is problematic and cannot be implemented 
correctly in the current framework. Rectifying this situation through a twisting of the 
symmetry seems problematic and a revision of the non-commutative relations themselves 
may be required. 

It will also be interesting to study, even in two dimensions, non-commutative 
fermion gases confined by a mutual attractive interaction between particles rather than 
a confining potential. In this case an interesting interplay between the excluded volume, 
thermal length and interaction length scales may be expected. 
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Our goal is to derive an approximate expression for tlie 2D g-potential 

da / de d{e, a) log[l + e'^^^-^^"))] (A.l) 
at low temperatures and high densities. Here /(a) = c + wa and 



(i(e, a) 



[a 



27re 



27re 



with 



e^{a) = a + 2 ± + 1. 
The values of a at which /(a) intersects the e±{Q.) curves are 



c — cw + 2w ±: 2-\/c — cw + 



(A.2) 
(A3) 

(A.4) 



(1-W)2 

See figure [H Let denote the e-integral in flA.ip for a fixed a. An expression for 
X(a) can be derived in terms of the three functions 



d{e' , Q;)de' = ed{e, a) + a;0(— a) H — arccos 

TT 



e_ + e+ — 2e 



e_ + e , - 4 



L(e, a) 



e'd{e', a)de 



2n 



IGtt 

/(e',a)de' = e/(e)-D(e) 



arccos 



arccos 



e+ - e_ 
e_ + e , + 2e - 4" 



e_ + e+ - 2e" 



(e_ + e+ - 2e) 



e(i(e, a) 



Here e± = e±{a). Important special cases are I{e-,a) = 0, /(e+,a) = 1 + a9(— a), 
L)(e_, a) = and a) = (e+ - e_)Vl6. 

To evaluate at low T (large b) there are three cases to consider: 

(i) /(a) > e+: Here — 6(e — /(a)) > and so exp[— 6(e — /(a))] is exponentially large 
in b. We approximate 

log[l + e-''(^-^("))] ^ -b{e - f{a)) (A.5) 

to produce 

X{a) ^ bfia)I{e+, a) - 6D(e+, a) = 6/(a)(l + a0(-a)) - ^b{e+ - e.f (A.6) 

In figure [U (a) this case corresponds to the region — 1 < a < 

(ii) e_ < f{a) < e+: Here the effective chemical potential /(«) lies inside the 
integration domain of the e integral and we use the Sommerfeld expansion [25] 
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to approximate the integral up to linear order in b ^. Integrating by parts twice 
and then expanding the integrand around e = f{a) produces 

X(a) ^ bLifia), a) + a). (A.7) 

In figure [T] (a) this case corresponds to the region a_ <«<«+. 

(iii) /(a) < e_: Here —b{e — f{a)) < and the integrand in exponentially small in b 
for all e G [e„,e+]. We therefore approximate ~ 0. In figured] (a) this case 
corresponds to the region a > «+. 

With X(a) known it remains to perform the a- integral in flA.ip . This involves 
extensive algebraic manipulations which we will not present here. The integrals 
themselves are however standard and may be found in texts such as [27]. 

References 

[1] Doplicher S, Fredenhagen K and Roberts J E 1995 Gommun. Math. Phys. 172 187 

[2] Douglas M R and Nekrasov N A 2001 Rev. Mod. Phys. 73 97 

[3] Scholtz F G, Gouba L, Hafver A and Rohwer C M 2009 Jnl. Phys. A 42 175303 

[4] Li K and Dulat S 2006 Eur. Phys. J. C 46 825 

[5] Vilela Mendesa R 2005 Eur. Phys. J. C 42 445 

[6] Bemfica F S and Girotti H O 2005 Jnl. Phys. A 38 L539 

[7] Khan S, Chakraborty B and Scholtz F G 2008 Phys. Rev. D 78 025024 

[8] Chaichian M, Sheikh- Jabbari M M and Tureanu A 2001 Phys. Rev. Lett. 86 2716 

[9] Chair N and Sheikh- Jabbari M M 2001 Phys. Lett. B 504 141 
[10] Liaoa Y and Dehneb C 2003 Eur. Phys. J. C 29 125 
[11] Ohl T and Renter J 2004 Phys. rev. D 70 076007 

[12] Garcfa-Compean H, Obregon O and Ramirez C 2002 Phys. Rev. Lett. 88 161301 
[13] Chaichian M, Kulish P P, Nishijima K and Tureanu A 2004 Phys. Lett. B 604 98 
[14] Alexander S, Brandenberger R and Magueijo J 2003 Phys. Rev. D 67 081301 
[15] Scholtz F G, Chakraborty B, Govaerts J and Vaidya S 2007 Jnl. Phys. A 40 14581 
[16] Wess J 2004 Deformed coordinate spaces: Derivatives arXiv:hep-th/0408080 
[17] Aschieri P, Blohmann C, Dimitrijevic M, Meyer F, Schupp P and Wess J 2005 Class. Quant. Grav. 
22 3511 

[18] Balachandran A P, Mangano G, Pinzul A and Vaidya S 2006 Int. J. Mod. Phys. A 21 3111 
[19] Chakraborty B, Gangopadhyay S, Hazra A G and Schohz F G 2006 Jnl. Phys. A 39 9557 
[20] Fiore G and Wess J 2007 Phys. Rev. D 75 105022 

[21] Fiore G 2008 On the consequences of twisted Poincare symmetry upon QFT on Moyal 

noncommutative spaces arXiv:hep-th/0809.4507 
[22] Schohz F G and Govaerts J 2008 Jnl. Phys. A 41 505003 
[23] Arriola E R and Dehesa J S 1989 Nuovo Gimento B 103 611 
[24] Dette H and Studden W J 1995 Gonstructive Approximation 11 227 
[25] Pathria R K 1972 Statistical Mechanics (Oxford: Pergamon) 
[26] Shariati A, Khorrami M and Fatollahi A H 2010 Jnl. Phys. A 43 285001 

[27] Gradshteyn I S and Ryzhik I M 2007 Table of integrals, series and products 7th ed (Burlington: 
Academic Press) 

[28] Chakraborty B, Kuznetsova Z and Toppan F 2010 J.Maths.Phys. 51 112102 

[29] Sinha D, Chakraborty B and Scholtz F G, "Non-commutative Quantum Mechanics in Three 
Dimensions and Rotational Symmetry", arXiv:hep-th/1108.2569 



